GAUSSIAN ESTIMATES FOR FUNDAMENTAL SOLUTIONS OF 
SECOND ORDER PARABOLIC SYSTEMS WITH 
TIME-INDEPENDENT COEFFICIENTS 



SEICK KIM 

Abstract. Auscher, Mcintosh and Tchamitchian studied the heat kernels 
of second order elliptic operators in divergence form with complex bounded 
measurable coefficients on R™. In particular, in the case when n = 2 they 
obtained Gaussian upper bound estimates for the heat kernel without imposing 
further assumption on the coefficients. We study the fundamental solutions of 
the systems of second order parabolic equations in the divergence form with 
bounded, measurable, time-independent coefficients, and extend their results 
to the systems of parabolic equations. 



1. Introduction 

In 1967, Aronson [I] proved Gaussian upper and lower bounds for the fundamen- 
tal solutions of parabolic equations in divergence form with bounded measurable 
coefficients. To establish the Gaussian lower bound Aronson made use of the Har- 
nack inequality for nonnegative solutions which was proved by Moser in 1964 (see 
1TJ ) . Related to Moser's parabolic Harnack inequality, we should mention Nash's 
earlier paper [18] where the Holder continuity of weak solutions to parabolic equa- 
tions in divergence form was established. In 1985, Fabes and Stroock [TO] showed 
that the idea of Nash could be used to establish a Gaussian upper and lower bound 
on the fundamental solution. They showed that actually such Gaussian estimates 
could be used to prove Moser's Harnack inequality. We note that Aronson also ob- 
tained Gaussian upper bound estimates of the fundamental solution without using 
Moser's Harnack inequality. 

In [2], Auscher proposed a new proof of Aronson's Gaussian upper bound esti- 
mates for the fundamental solution of second order parabolic equations with time- 
independent coefficients. His method relies crucially on the assumption that the 
coefficients are time-independent and thus it does not exactly reproduce Aronson's 
result, which is valid even for the time-dependent coefficients case. However, his 
method is interesting in the sense that it carries over to equations with complex 
coefficients provided that the complex coefficients are a small perturbation of real 
coefficients. Along with this direction, Auscher, Mcintosh and Tchamitchian also 
showed that the heat kernel of second order elliptic operators in divergence form 
with complex bounded measurable coefficients in the two dimensional space has a 
Gaussian upper bound (see [3] and also [5]). 

We would like to point out that a parabolic equation with complex coefficients 
is, in fact, a special case of a system of parabolic equations. From this point of 
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view, Hofmann and the author showed that the fundamental solution of a parabolic 
system has an upper Gaussian bound if the system is a small perturbation of a di- 
agonal system, which, in particular, generalized the result of Auscher mentioned 
above to the time-dependent coefficients case (see [H]). However, the above men- 
tioned result of Auscher, Mcintosh and Tchamitchian regarding the heat kernel of 
two dimensional elliptic operators with complex coefficients does not follow directly 
from our result. 

One of the main goals of this article is to provide a proof that weak solutions of 
the parabolic system of divergence type with time-independent coefficients associ- 
ated to an elliptic system in two dimensions enjoy the parabolic local boundedness 
property and to show that its fundamental solution has a Gaussian upper bound. 
More generally, we show that if weak solutions of an elliptic system satisfy Holder 
estimates at every scale, then weak solutions of the corresponding parabolic system 
with time-independent coefficients also satisfies similar parabolic Holder estimates 
from which, in particular, the parabolic local boundedness property follows easily. 
Also, such an argument allows one to derive Holder continuity estimates for weak 
solutions of parabolic equations with time-independent coefficients directly from 
De Giorgi's theorem [7] on elliptic equations, bypassing Moser's parabolic Harnack 
inequality. In fact, this is what Auscher really proved in the setting of complex 
coefficients equations by using a functional calculus method (see [5] and also [I], 
[5]). Even in those complex coefficients settings, we believe that our approach is 
much more straightforward and thus appeals to wider readership. 

Finally, we would like to point out that in this article, we are mainly interested 
in global estimates and that we do not attempt to treat, for example, the systems 
with lower order terms, etc. However, let us also mention that, with some extra 
technical details, our methods carry over to those cases as well as to the systems 
of higher order; see e.g. [4], [5] for the details, and also Remark 13.51 

The remaining sections are organized in the following way. In Section [2] we give 
notations, definitions, and some known facts. We state the main results in Section[3] 
and give the proofs in Section |U 

2. Notation and definitions 

2.1. Geometric notation. 

(1) 1" = n-dimensional real Euclidean space. 

(2) x — (x\, ■ ■ ■ ,x n ) is an arbitrary point of R™. 

(3) X — (x, t) denotes an arbitrary point in R n+1 , where x £ R™ and tel. 

(4) B r (x) — {y G l n : \y — x\ < t} is an open ball in R™ with center x and 
radius r > 0. We sometimes drop the reference point a; and write B r for 
B r (x) if there is no danger of confusion. 

(5) Q r {X) = {(y, s) e R" +1 : \y - x\ < r and t - r 2 < s < t}. We sometimes 
drop the reference point X and write Q r for Q r (X). 

(6) Q*(X) = {(y,s) e R n+1 : \y - x\ < r and t < s <t + r 2 }. 

(7) Qr, s { X ) = {(y,s) eQ r (X)}- i.e., Q T ,,{X) = B r (x) x {s} if s e (t-r 2 ,t) 
and Q riS (X) — otherwise. We sometimes drop the reference point X and 
write Q r ,s for Q r ^ s (X). 

(8) For a cylinder Q — £1 x (a, b) C R™ +1 , dpQ denotes its parabolic boundary, 
namely, dpQ = dtt x (a.b) U il x {a}, where 9f2 is the usual topological 
boundary of £1 C R™ and is its closure. 
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2.2. Notation for functions and their derivatives. 

(1) For a mapping from Q C M" to R N , we write f(x) = {^(x), f N (x)) T 
as a column vector. 

(2) /q = j^y Jq /, where |Q| denotes the volume of Q. 

(3) u t = du/dt. 

(4) D^u = DiU = u Xi = du/dxi. 

(5) Du — (u Xl , . . . , u Xn ) T is the spatial gradient of u = u(x, t). 

(6) For / = (,f \ . . . , f N ) T , Df = (Df\ . . . , Df N ); that is Df is the n x N 
matrix whose i-th column is Df 1 . 

2.3. Function spaces. 

(1) For £1 C K" and p > 1, L p (f2) denotes the space of functions with the 
following norms: 

\\u\\ LP{n) = (^jju(x)\ p dx^j and ||u|| ioo(n) = esssup a |m| . 

(2) C M (f2) denotes the space of functions that are Holder continuous with the 
exponent /x E (0, 1], and 

Mc*.(n) = sup < oo. 

(3) The Morrey space M 2,M (r2) is the set of all functions u E L 2 (f2) such that 

1/2 

_ ,,„„ f / L.|2 ^ 

lilf 2 if(Si) 



|2 

B P (i)C!) V Jb„(x) 



W u \\ M ^(m = sup p " / |«r I < oo. 



(4) Cp(Q) denotes the space of functions defined on Q C R™ +1 such that 
Mc£(Q) = sup - < oo, 

where d P (X,X') = max (|a; - x'\ , ^/|t - t'\j. 
2.4. Elliptic and parabolic systems and their adjoints. 

Definition 2.1. We say that the coefficients A"? (x) satisfy the uniform ellipticity 
condition if there exist numbers vq, M > such that for all x E R" we have 

a P(„\t. c \ -> It I 2 / 



(2.1) {A ap {x)^ a )>u \^ and ^(z)^, r/ Q ) < M \i\ \ V \ , 

where we used the following notation. 

(1) For a, (3 = 1, . . . ,n, A af3 (x) are N x N matrices with (i, j)-entries A°f(x). 

n n „ 

(2) €„ = (^,---,^) r and |£| 2 = £ £ |£| . 

a— 1 i= 1 

(3) (A at3 {x)i ,r)o) = E E AfW&i. 

We emphasize that we do not assume that the coefficients are symmetric. 
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Definition 2.2. We say that a system of N equations on M. n 

N n 

E E ^(<w^" j ) = o (i = i,...,N) 

3=1 o,/3=l 

is elliptic if the coefficients satisfy the uniform ellipticity condition. We often write 
the above system in a vector form 

n 

(2.2) Lu := E A*(^ a/S (aO-D/J«) = 0, u = {u 1 . . . ,u N ) T . 

a, 13=1 

The adjoint system of (|2.2|) is given by 

n 

(2.3) L*u := E A, ((A°^)*(a;)D fl ti) - 0, 

a, j3=l 

where (A Q/3 )* = (A' 3 ") 7 , the transpose of A 13 ". 

Definition 2.3. We say that a system of N equations on R ra+1 

N n 

«!"EE D Xa {Af{x)D x y)=Q (i = l,...,N) 

j=l a,(3 = l 

is parabolic if the (time-independent) coefficients satisfy the uniform ellipticity 
condition. We often write the above system in a vector form 

n 

(2.4) u t -Lu :=u t - E D a (A a0 (x)D u) = 0. 

a,P=l 

The adjoint system of (|2.4[) is given by 

n 

(2.5) ut + i*« := u t + E D « (O 4 ^)*^)- /^) = °! 

a,/3=l 

where (A Q/3 )* = (A' 3 ") 7 , the transpose of A Pa . 

2.5. Weak solutions. In this article, the term "weak solution" is used in a rather 
abusive way. To avoid unnecessary technicalities, we may assume that all the coef- 
ficients involved are smooth so that all weak solutions are indeed classical solutions. 
However, this extra smoothness assumption will not be used quantitatively in our 
estimates. This is why we shall make clear the dependence of constants. 

(1) We say that it is a weak solution of (|2.2p in £1 C R™ if it is a (classical) 
solution of (|2.2|) in SI and u, Du E L 2 (fl). 

(2) We say that u is a weak solution of (|2.4|) in a cylinder Q = Qx (a, b) C M™ +1 
if u is a (classical) solution of (|2.2[) in Q and it, Du S L 2 (Q), u(-, t) s £ 2 (il) 
for all a < £ < &, and sup Q<t<{) i)||i2(n) < °°- 

2.6. Fundamental solution. By a fundamental solution (or fundamental matrix) 
T(x, t ; y) of the parabolic system (|2.4p we mean an iV x TV matrix of functions defined 
for t > which, as a function of (a;, t), is a solution of (12. 4[) (i.e., each column is a 
solution of l|2.4p ). and is such that 

(2.6) lim f T(x,t;y)f(y)dy = f(x) 
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for any bounded continuous function / = (/ , . . . ,f N ) T , where T(x,t;y)f(y) de- 
notes the usual matrix multiplication. 

2.7. Notation for estimates. We employ the letter C to denote a universal con- 
stant usually depending on the dimension and ellipticity constants. It should be 
understood that C may vary from line to line. We sometimes write C — C(a, (3, . . .) 
to emphasize the dependence on the prescribed quantities a, (3, . . .. 

2.8. Some preliminary results and known facts. 

Lemma 2.4 (Energy estimates). Letu be a weak solution of (|2.4p inQ R = Q R {X). 
Then for < r < R, we have 

sup / \u{-,s)\ 2 +f \Duf < ° [ \u\ 2 . 

Proof. See e.g., [H Lemma 2.1, p. 139]. □ 

Lemma 2.5 (Parabolic Poincare inequality). Let u be a weak solution of (|2.4| in 
Qr = Qr{X). Then there is some constant C — C(h,Mq) such that 

\u-u Qn \ 2 < CR 2 I \Du\ 2 . 

JQr 

Proof. See e.g., PjJJ Lemma 3]. □ 

Lemma 2.6. Let Q2R — Q2r{Xq) be a cylinder in K n+1 . Suppose u £ L 2 (Q2r) 
and there are positive constants fj, < 1 and M such that for any X S Qr and any 
r G (0, R) we have 

I \u-u Qr{x) \ 2 <M 2 r n+2+2 ». 

JQr(X) 

Then u is Holder continuous in Qr with the exponent /j, and Mc^(Q„) ^ C( n i (J,)M. 

Proof. See e.g., [HI Lemma 4.3, p. 50]. □ 

Definition 2.7 (Local boundedness property). We say that the system (|2.4p sat- 
isfies the local boundedness property for weak solutions if there is a constant M 
such that all weak solutions u of (|2.4j) in Q2r(X) satisfy the estimates 

sup |u| < M (-— / \u\ 2 . 

Qr(X) \\Q2r\ JQ 2T (X) J 

Similarly, we say that the adjoint system (|2.5p satisfies the local boundedness prop- 
erty if the corresponding estimates hold for weak solutions u of (|2.5|) in Q2 r (X). 

Theorem 2.8 (Theorem 1.1, [12]). Assume that the system (|2.4[) and its adjoint 
system (|2.5|) satisfy the local boundedness property for weak solutions. Then the 
fundamental solution of the system (|2.4I) has an upper bound 



(2-7) |r(x,t;y)| op <C t-"/ 2 exp(-^ 



where \T(x, t] y)\ denotes the operator norm of the fundamental matrix T(x, i; y). 
Here, Cq = Co(n, Vq, Mo, M) and kg = ko(fO) Ma)- 
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3. Main results 



Definition 3.1. We say that an elliptic system (J272J) satisfies the Holder estimates 
for weak solutions at every scale if there exist constants fj,o > and Hq such that 
all weak solutions u of the system in Bi r = Bi r {xo) satisfy the following estimates 

(3-1) Mc7«( Br ) < H r-W^ \\u\\ L2{B2r) . 

Similarly, we say that a parabolic system (12. 4|) satisfies Holder estimates for weak 
solutions at every scale if there exist constants [i\ > and H± such that all weak 
solutions u of the system in Qi r = Qir{Xo) satisfy the following estimates 

(3-2) Mc« (Qr ) < H ir -W™ + rt \\u\\ L2(Q2r) . 

Remark 3.2. Elliptic systems with constant coefficients satisfy the above property, 
and in that case, the ellipticity condition (|2.ip can be weakened and replaced by the 
Legendre-Hadamard condition. De Giorgi's theorem [7] states that the property is 
satisfied if N = 1. The property is also satisfied if n = 2 and it is due to Morrey 
(see Corollary 13. 6[) . Some other examples include, for instance, a certain three 
dimensional elliptic system which was studied by Kang and the author in [13) . 

We shall prove the following main results in this paper: 

Theorem 3.3. If an elliptic system (|2.2p satisfies the Holder estimates for weak 
solutions at every scale, then the corresponding parabolic system (|2.4|) with time- 
independent coefficients also satisfies the Holder estimates for weak solutions at 
every scale. 

Theorem 3.4. Suppose that the elliptic system (|2.2p and its adjoint system (|2.3|) 
defined on K™ both satisfy the Holder estimates for weak solutions at every scale 
with constants /j,q,Hq. Let T(x.t;y) be the fundamental solution of the parabolic 
system (|2.4p with the time-independent coefficients associated to the elliptic system 
(|2.2p . Then T(x, t;y) has an upper bound 




(3.3) |r(;M;y)| op <Coi-™/ 2 exp^-^: 

where Cq — Co(n, vq, Mq, //q, Ho) and kg — ko(i/Q, Mq). Here, \T(x, t; y)\ op denotes 
the operator norm of fundamental matrix T(x,t;y). 

Remark 3.5. We would like to point out that (|3 . 3j) is a global estimate. Especially, 
the bound (|3.3p holds for all time t > 0. Suppose that the elliptic system (|2.2p and 
its adjoint system (|2.3p enjoy the Holder estimates for weak solutions up to a fixed 
scale i?o; that is, there is a number Rq > such that if u is a weak solution of 
either (|2.2p or (|2.3p in B r = B r (x) with < r < i?o, then u is Holder continuous 
and satisfies 

[«lc«( Bf) <H r-^ 2 +^\\u\\ L2{B2r) . 
Then, the statement regarding the bound (|3.3|) for the fundamental solution should 
be localized as follows: For any given T > 0, there are constants ko = k {v ,M ) 
and C = C (n, v , M , fx ,H ,Ro, T) such that (O holds for < t < T. 

Corollary 3.6. Let T(x, t; y) be the fundamental solution of the parabolic system 
()2.4p with time-independent coefficients associated to an elliptic system (|2.2p defined 
on R 2 . Then T(x, t; y) has an upper bound p.3p with the constants Co, ko depending 
only on the ellipticity constants Uq,M. 
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Proof. First, let us recall the well known theorem of Morrey which states that any 
two dimensional elliptic system (|2.2|) with bounded measurable coefficients satisfies 
the Holder estimates for weak solutions at every scale, with the constants /j,o,Hq 
depending only on the ellipticity constants (see, [16l pp. 143-148]). Next, note that 
the ellipticity constants Vq,Mq in ()2.1|) remain unchanged for A ff (x ) — Aj™(x). 
Therefore, the corollary is an immediate consequence of Theorem 13.41 □ 

Remark 3.7. In fact, the converse of Theorem l2.8l is also true (see [12l Theorem 1.2]). 
Therefore, in order to extend the above corollary to the parabolic system with 
time-dependent coefficients, one needs to show that the system satisfies the local 
boundedness property for weak solutions. Unfortunately, we do not know whether 
it is true or not if the coefficients are allowed to depend on the time variable. If 
n > 3, it is not true in general, even for the time-independent coefficients case since 
there is a famous counter-example due to De Giorgi (see [8]). 



4. Proof of Main Results 
4.1. Some technical lemmas and proofs. 

Lemma 4.1. If u is a weak solution of the parabolic system with time-independent 
coefficients (|2 .4[) in Qr = Qr(Xq), then u t G L 2 {Q r ) for r < R and satisfies the 
estimates 

(4-1) \\u t \\ L2(Qr) <C{R-r)-'\\Du\\ L2(QR) . 

In particular, if u is a weak solution of (|2.4[) in Qi r , then the above estimates 
together with the energy estimates yield 

(4-2) WM LHQr) <Cr- 2 \\u\\ L2(Q2r) . 

Proof. We first note that if the coefficients are symmetric, (i.e., A"? = Aj") this 
is a well known result; a proof for such a case is found, for example, in [14j pp. 
172-181] or in [9l pp. 360-364]. However, the standard proof does not carry over to 
the non-symmetric coefficients case and for that reason, we provide a self-contained 
proof here. 

Fix positive numbers a, r such that a < r < R. Let ( be a smooth cut-off 
function such that £ = 1 in Q a , vanishes near dpQ T , and satisfies 

< C < 1 and \( t \ + \D(\ 2 <C{r~a)- 2 . 

Note that on each slice Q T , S , we have 

= ^ (u t -D a (A al3 D p u)) -C 2 u t 
= [ C 2 \u t \ 2 + I ( 2 (A al3 D l3 u,D a u t )+ [ 2((A a(3 D p u,D a (u t ). 



8 



SEICK KIM 



Therefore, we find by using the Cauchy-Schwarz inequality that 

/ Q 2 \u t \ 2 <c( C, 2 \Du\\Du t \ + C ( (\Du\\D(\\u t 

JQt . J Qt , J Qt . 



<-/ C 2 \Du t \ 2 + -[ C 2 \Du\ 2 + cf \DC\ 2 \Du\ 2 

>Qt,s JQr.s 



2 

Thus we have 



2.1,, 
1 



C 2 M 



(4.3) f C 2 \u t \ 2 <J C 2 \Du t \ 2 +Cf C 2\ Du f + c [ \ DC \ 2 \Du\ 2 . 

J q t JQr 6 Jq t JQt 

Since u t also satisfies (|2.4|) . the energy estimates yield 

(4-4) / C 2 \Du t \ 2 <--^ ( \u t \ 2 . 

JQt ( T - ) *> Qt 

This is the part where we exploit the assumption that the coefficients are time- 
independent. Combining (|4.3p and (|4.4|) . we have 



\u t \ 2 <-^- 2 [ |n t | 2 + -/ \Du\* + T °[ \Du\ 2 . 

'Q a ( T - a) JQt e JQt ( t - ) JQt 

If we set e = (r — cr) 2 /2Co, we finally obtain 

k|2 ^/ ,"" |2+ f^fX. |d "^ 

Here, we emphasize that C is a constant independent of <r, r. Then by a standard 
iteration argument (see e.g. [Ill Lemma 3.1, pp. 161] ) , we have 

(4.5) / |tt t | 2 < , D C .- / \Du\ 2 forO<r<i?. 

The proof is complete. □ 

Lemma 4.2. If u is a weak solution of the parabolic system with time-independent 
coefficients (12. 4p in Q2r = Qir(Xo), then Du(-, s),u t (-, s) G L 2 (Q rs ) for all s £ 
[to — r 2 , to], and satisfy the following estimates uniformly in s S [to — r 2 , 

(4-6) ||£>tx(-, S )|| L2(0ris) <Cr- 2 ||«|| L2(Q2r) , 

(4-7) \\u t (;s)\\ L2{Qrs) <Cr- 3 \\u\\ L2(Q2r) . 

Proof. By the energy estimates applied to Ut we obtain 

(4.8) sup / \u t (-, s)\ 2 < £L f \u t \ 2 . 

to-r 2 <s<t jQ r , s r JQ 3r/2 

On the other hand, the estimates (|4. 5[) and the energy estimates (this time, applied 
to u itself) yield 

(4.9) / \u t \ 2 <^[ \Du\ 2 <-( \u\ 2 . 
Combining (|4.8[) and (|4.9[) together, we have the estimates (|4.7p . 



GAUSSIAN ESTIMATES 



9 



Next, assume that u is a weak solution of (|2.4|) in Q^ r = Q4 r (Xo). Let £ be a 
smooth cut-off function such that £ = 1 in Q r , vanishes near dpQ2r, and satisfies 

(4.10) < C < 1 and | C* I + |£>C| 2 < Cr~ 2 . 

Note that on each slice Q2r,s, we have 

0= / (u t ~D a (A afi D u)) -C 2 u 

JQ2r,s ^ ' 

= J ( 2 u t -u + J ( 2 (A af3 D f3 u,D a u) + 2c(A afi D f3 u,D a (u 
Using the ellipticity condition and the Cauchy-Schwarz inequality, we find 



f ( 2 \Du\ 2 <f ( 2 \u t \\u\+C [ (\Du\\D(\\u\ 



V ° f ( 2 \Du\ 2 . 



2 

Then by (I4.10p . (|4.7|) . and the energy estimates, for all s £ [(q - r 2 ,t ], we have 
\Du\ 2 <e[ |« t | 2 + -/ H 2 + ^f \u\ 2 

<— / l«l +— / H+-/ l«l • 

If we set e = r 2 , then the above estimates (|4. 1 1[) now become 

\Du\ 2 <-( \u\ 2 , 

r JQir 

from which the estimates (I4.6[) follows by a well known covering argument. □ 

Lemma 4.3. Assume that the elliptic system <\2.2\ satisfies the Holder estimates 
for weak solutions at every scale with constants /iq,Hq. Let u be a weak solution of 
the inhomogeneous elliptic system 

(4.12) D a (A aP (x)D p v) = f in B 2 = B 2 (x ), 

where f belongs to the Morrey space M 2 ' X (B 2 ) with A > 0. Then, for any 7 > 
with 7 < 70 = min(A + 4, n + 2/iq) (we may take 7 = 70 if 70 < n) there exists 
C = C(n, vq, Mq, ho, Hq, X, 7) such that u satisfies the following local estimates 



(4.13) 



/ \Du\ 2 <C (r"<- 2 f \Du\ 2 +r"'- 2 \\f\ 
JB T (x) \ Jb 2 



M 2 - X (B 2 ) 



uniformly for all x S B\ = Bi(xq) and < r < 1. Moreover, if 7 < n, then u 
belongs to the Morrey space M 2,1 (Bi) and 

(4.14) \\u\\ M ^ [Bl) < C (||u|| £a(Ba) + \\Du\\ L2[B2) + ||/|| M2 ,x (S2) ) . 
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Proof. First, we note that the property (|3.1|) implies that for all < p < r and 
x G R™, we have 



L 



\Du\ 2 <C-H (B.) / \Du\ 2 



'B p (x) xr/ JB T {x) 

In the light of the above observation, the estimates (|4. 13|) is quite standard and is 
found, for example, in [111 Chapter 3]. Then, by Poincare inequality we have 



2 

M 2 ^{B 2 ) 



(4-15) f \u-u Br(x) \ 2 <Cr^ (\\Du\\ 2 L2{B2) + \\f\ 

JB r (x) V 

uniformly for all x G B\ = -Bi(O) and < r < 1. It is well known that if 7 < n, 
then the estimates (|4.15l) yield (|4. 14|) (see e.g. pU Chapter 3]). □ 



4.2. Proof of Theorem 13.31 Let it be a weak solution of (|2.4|) in a cylinder 
Q4 = Qa(0)- We rewrite (|2 .4[) as Lu = u t . By Lemma T4. 21 we find that u t (-,s) is 
in L 2 {Q2 : s) and satisfies 

W u t(-,s)\\ L 2 {Q2 s) < C\\u\\ L2(Qi) for all - 4 < s < 0. 

Therefore, we may apply Lemma 14.31 with f = Ut and A = 0, and then apply 
Lemma |4~21 to find that for all x G £?i(0) and < r < 1, we have 

/ S )| 2 < Cr^ 2 S )|| 2 2(Q2 + || Mt (, S )|| 2 2(Q2 } ) 

(4.16) Jb t {x) v ' v ' ' ' v ' 

< Cr~<- 2 \\u\\ 2 L2{Qi) uniformly in s G [-4, 0] 

for all 7 < min(4, n + 2fio). 

By Lemma [2751 and then by (|4. 16|) we find that for all X = (x, t) G Q\ and r < 1 

/ |m -Mq,.(x)| 2 < C*r 2 /" / \Du(y,s)\ 2 dyds 

(4.17) JQr(X) Jt-r 2 JB r {x) 

<Cr 2 +-'\\u\\ 2 L2(Qi) . 

Note that if n < 3, then we may write 7 = n + 2/i for some > 0. In that case, 
(|4.17p now reads 

(4-18) / \u-u QAX) \ 2 < CV»+ 2 + 2 " ||n||i 2(Q4) 



for all X G Qi and r < 1. Therefore, if n < 3, then Lemma 12761 yields the estimates 
(4-19) ^H { Q 1/2 )<C\\u\\ L2(Qi) . 

We have thus shown that in the case when n < 3, any weak solution u of (|2.4| in a 
cylinder Q 4 = (54(0) satisfies the above a priori estimates (|4.19[) provided that the 
associated elliptic system satisfies the Holder estimates for weak solutions at every 
scale. The general case is recovered as follows. For given X = (x ,t Q ) and r > 0, 
let us consider the new system 

n 

(4.20) u t -Lu:=u t - ^ D a (A af} {x)D l3 u) = 0, 

u,f3=l 

where A a P(x) — A a/3 (xo + rx). Note that the associated elliptic system Lu = 
also satisfies the Holder estimates for weak solutions at every scale. Moreover, the 
ellipticity constants vq, Mq remain the same for the new coefficients A a @ . Let u be 
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a weak solution of (12. 41) in Q4 r (Xo). Then u(X) = u(x, t) :— u(xq + rx, to + r 2 t) 
is a weak solution of (|4.20|) in Q4(0) and thus u satisfies the estimates (|4.19j) . By 
rescaling back to Q4 r (Xo), the estimates (I4.19|) become 

(4-21) M^(Q. /2 )<^- (,l/2+1+Al) ll«ll^(Q 4 .)- 

Thus, when n < 3, the theorem now follows from a well known covering argument. 

In the case when n > 4, we invoke a bootstrap argument. For the sake of 
simplicity, let us momentarily assume that 4 < n < 7. Let u be a weak solution 
of (|2.4p in Qs = Qs(0)- Let us fix Xq = (xo,to) € Q2(0) and observe that Ut also 
satisfies the system (12.41) in Q^Xq). Thus, by a similar argument that led to (|4.16l) . 
we find that for all x <E Bi(xo) and < r < 1 we have 

(4.22) f \Du t (;.s)\ 2 < Cr'- 2 \\u t \\ 2 L2(Qi(Xn)) uniformly in s e [to -4, to], 

for all 7 < 4 (we may take 7 = 4 if n > 4). Then, by (|4. 14(1 in Lemma |4.3[ 
Lemma FTTl and Lemma T4. 2 1 we conclude that 

(4.23) ll«t(-.s)|| M2 , 7(Bl(:Bo)) < C\\u\\ L2{Qsm for all s€ [t -4,t ]. 

Since the above estimates (|4.23|) hold for all Xq = (xo, to) € Q2(0), we find that, in 
particular, u t (-,s) belongs to M 2 ' 7 (i?2(0)) for all —4 < s < 0, and satisfies 

(4.24) \\M;s)\\ M 2, HB2m <C\\u\\ L2{Qa(0)) for all se [-4,0], 

where we also used (14. 7[) of Lemma 14.21 

The above estimates (|4.24p for u t now allows us to invoke Lemma l4T31 with f = u t 
and A — 7. Then, by Lemma [4.31 and Lemma [4.2[ we find that for all x E Si(0) 
and < r < 1, we have 

[ \Du(;.s)\ 2 < Cr^- 2 \\u\\ 2 L2{Qsm uniformly in s e [-4,0] 

J B r (x) 

for all 7 < min(7+4, n+2/io). Since we assume that n < 7, we may write 7 = n+2/Z 
for some 7i > 0. By the exactly same argument we used in the case when n < 3, 
we derive the estimates 

Mc£(Q 1/2 ) <C\\u\\ L2{Qs) , 

and the theorem follows as before. 

Finally, if n > 8, we repeat the above process; if u is a weak solution of (|2.4p in 
Qi6(0), then u t (-, s) is in M 2,7 (_B 1 (0)) for all 7 < 8 and so on. The process cannot 
go on indefinitely and it stops in k = [n/4] + 1 steps. The proof is complete. □ 

4.3. Proof of Theorem 13.41 The proof is based on Theorem 12.81 the proof of 
which, in turn, is found in [12] . By Theorem 12.81 we only need to establish the 
local boundedness property for weak solutions of the parabolic system (12.4[) and for 
those of its adjoint system (|2.5p . 

From the hypothesis that the elliptic system (|2.2p satisfies the Holder estimates 
for weak solutions at every scale, we find, by Theorem 13.31 that the parabolic 
system (12.41) with the associated time-independent coefficients also satisfies the 
Holder estimates for weak solutions at every scale; that is, there exist some constants 
yU > and C, depending on the prescribed quantities, such that if u is a weak 
solution of (12.41) in Q^ r (X), then it satisfies the estimates 

M^(Q 2 ,.)<c^ ( " /2+1+M) ll«ll i2(Q40 - 
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Let us fix Y £ Q r = Qr(X). Then, for all Z £ Q r {Y) C Q 2 r(X), we have 

(4.25) \u(Y)\ < \u(Z)\+d P (YZr-[u] c , {Q2r) < \u(Z)\+Cr-^ 2+ ^ \\u\\ LHQir) . 
By averaging (|4.25p over Q r [Y) with respect to Z, we derive (note \Q r \ — Cr n+2 ) 

\u(Y)\ < Cr"(™+ 2 ) ||«|| L1(Q ,. (y)) + Cr-(»/2+D \\ u \\ L2{Qir} . 

Since Y £ Q r is arbitrary, we find, by Holder's inequality, that u satisfies 

HI L - (<W < Cr-^+V \\ u \\ LHQir) 

for some constant C — C(n, Mo, fj-o, Hq). 

To finish the proof, we also need to show that if u is a weak solution of the adjoint 
system (|2.5[) in Q\ r — Q\ r {X), then it satisfies the local boundedness property 

(4.26) ||«|| LOO(Q , ) <Cr~("/ 2 + 1 )|| U || i2(Q , r) . 

The verification of (I4.26|) requires only a slight modification of the previous argu- 
ments (mostly, one needs to replace Q r by Q* and so on), but it is rather routine 
and we skip the details. □ 



References 

1. Aronson, D. G. Bounds for the fundamental solution of a parabolic equation. Bull. Amer. 
Math. Soc. 73 (1967) 890-896. 

2. Auscher, P. Regularity theorems and heat kernel for elliptic operators. J. London Math. Soc. 
(2) 54 (1996), no. 2, 284-296. 

3. Auscher, P.; Mcintosh, A.; Tchamitchian, Ph. Heat kernels of second order complex elliptic 
operators and applications. J. Funct. Anal. 152 (1998), no. 1, 22-73. 

4. Auscher, P.; Qafsaoui, M. Equivalence between regularity theorems and heat kernel estimates 
for higher order elliptic operators and systems under divergence form. J. Funct. Anal. 177 
(2000), no. 2, 310-364. 

5. Auscher, P.; Tchamitchian, Ph. Square root problem for divergence operators and related 
topics. Asterisque No. 249 (1998) 

6. Davies, E. B. Explicit constants for Gaussian upper bounds on heat kernels, Amer. J. Math. 
109 (1987), no. 2, 319-333. 

7. De GiORGI, E. Sulla differenziabilita e V analiticita delle estremali degli integrali multipli 
regolari. (Italian) Mem. Accad. Sci. Torino. CI. Sci. Fis. Mat. Nat. (3) 3 (1957), 25-43. 

8. De GiORGI, E. Un esempio di estremali discontinue per un problema variazionale di tipo 
ellittico. (Italian) Boll. Un. Mat. Ital. (4) 1 (1968), 135-137. 

9. Evans, L. C. Partial differential equations. American Mathematical Society, Providence, RI, 
1998. 

10. Fabes, E. B.; Stroock, D. W. A new proof of Moser's parabolic Harnack inequality using the 
old ideas of Nash. Arch. Rational Mech. Anal. 96 (1986), no. 4, 327-338. 

11. Giaquinta, M. Multiple integrals in the calculus of variations and nonlinear elliptic systems. 
Princeton University Press:Princeton, NJ, 1983. 

12. Hofmann, S.; Kim, S. Gaussian estimates for fundamental solutions to certain parabolic 
systems. Publ. Mat. 48 (2004), 481-496. 

13. Kang, K.; Kim, S. On the Holder continuity of solutions of a certain system related to 
Maxwell's equations. SIAM J. Math. Anal. 34 (2002), no. 1, 87-100 (electronic). 

14. Ladyzenskaja, O. A.; Solonnikov, V. A.; Ural'ceva, N. N. Linear and quasilinear equations of 
parabolic type. American Mathematical Society: Providence, RI, 1967. 

15. Lieberman, G. M. Second order parabolic differential equations. World Scientific Publishing 
Co., Inc., River Edge, NJ, 1996 

16. Morrey, C. B., Jr. Multiple integrals in the calculus of variations. Springer- Verlag New York, 
Inc., New York, 1966 

17. Moser, J. A Harnack inequality for parabolic differential equations. Comm. Pure Appl. Math. 
17 (1964) 101-134; Correction: Comm. Pure Appl. Math. 20 (1967) 231-236. 



GAUSSIAN ESTIMATES 



13 



18. Nash, J. Continuity of solutions of parabolic and elliptic equations. Amer. J. Math. 80 (1958) 
931-954. 

19. Struwe, M. On the Holder continuity of bounded weak solutions of quasilinear parabolic sys- 
tems. Manuscripta Math. 35 (1981), no. 1-2, 125-145. 

Mathematics Department, University of Missouri, Columbia, Missouri 65211 
E-mail address: seickamath.missouri.edu 



